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1.0 Introducticn 



This report describes sane preliminary results of an effort to relate 
total aerosol optical depth at visible wavelengths to the visibility in the 
marine boindary la^er (MBL) . Wfe present an atmospheric model which explains 
the observed height dependence of the aerosol extinction, examine aircraft 
and ship data from the MAGAT experiment, illustrate the influence of synop- 
tic scale weather patterns on aerosol structure and present a simple model 
to relate to total vertical depth. The ultimate goal of the project is 
to allow an estimation of IR ranges (due primarily to atmospheric water 
vapor density, Q) from satellite remote sensing data. A feasibility study 
on the concept was performed on two NFS data sets (JASIN and CEWCQ^ 78) us- 
ing the relationship 

Q/Qg = a + b exp [-1.25 (0) 

where is the saturation water vapor density at 13°C, the surface 
extinction coefficient in and U is the wind speed in m/s. A least 

squares fit for 834 points yielded a = 1.02, b = -.27 and a correlation of 
0.56. 
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2. 0 Background 

2. 1 Optical Depth 

Let the atmo^heric trananissicn coefficient, T, due solely to 
aerosols be represented by 

T = exp (-U) (1) 

vt^iere 



V 



00 




a(Z)dZ 



( 2 ) 



o 

and a is the aerosol extincticn coefficient, Z the vertical coordinate, and 



a = j dN/dr Q(^,n)dr (3) 

o 

with r the aerosol particle radius, dN/dr the aerosol nunber density spec- 
trun and Q(^,n) the Mie scattering efficiency. If dN/dr is ]<nov\n frctn the 
surfece to the top of the atmosphere, then w can be calculated. Alter- 
nately, u can be measured fran satellites. The question is how is u re- 
lated to the surface visibility (extinction) . 

(insider an idealization of the atmospheric extinction coefficient 
vertical profile (Big. 1) 

V = <a>^ h + (4a) 

L (4b) 

vihere is the average extinction in a vvell mixed boundary la^er of 

depth, h, is the extinction in the clearer air just above the boundary 
la^er, L is an effective potential scale length for the decay of “( Z) for 
Z > h, is the optical depth of the troposphere above the mixed layer and 
is the optical depth of the stratosphere. This idealization is a result 



7 




Extinction 



Figure 1. Schematic representation of the height dependence of aerosol 

extinction coefficient, a, within and above the marine boundary 
layer of depth h. 
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of the typical veil-mixed nature of the marine boundary layer (Section 

2.2). The surfece extinction, “ , is related to <“> by the fector f 

o m 

“o = <5) 

vhich is less than cne due to the vertical structure of temperature and 

vater vapor in the marine boundary layer. 

If all of the relevant variables can be estimated, then ve can calcu- 
late using 

= f(u - \ - Wg)/h (6) 

Ebr example, and could be obtained fron the LOWTRAN Tropospheric 
Model , h frcm synoptic climatology and f established by model calculation 
and snpirical measurement. Typical values of L are 2 to 5 kn, might be 
10~^ to 10~^ in the visible. 

2.2 NFS Aerosol Boundary layer Model Structure 

The marine boundary la^er meteorological structure is illustrated in 
Fig. 2. The aerosol structure in the veil-mixed ccxitext is depicted in 
Fig. 3. The total aerosol volune spectral density, V(r^) , is the sun of 
continental, V^, and sea-salt, V^, conpcnents, where (V(r^) = dv/dr^) 



V = V+ V Z<h 

s c 

V = V Z > h 

c 

at sane reference hunidity (ve lave assuned no clouds) . 
is related to the nunber density by 
V(r) = 4/3 ifr^ dN/dr 



(7a) 

(7b) 

The volune density 
( 8 ) 



and the particle radius at saturation S is given by 

r = r^ g(S) (9) 

Tie specification of and is based cn onpirical parameteriza- 
tions. For example, can be represented as a Junge distribution and 
as a modified gaiima function, F (Goroch et al, 1982) 
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Figure 2. Marine boundary layer extinction for virtual potential tempera- 
ture, 9y, and water vapor density, Q, in the mixed-layer model 
structure. 




0 

Aerosol Volume Mixing Ratio 
at Reference Relative Humidity 



Figure 3. Similar to Fig. 2 but for the aerosol density. 
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(10a) 




Vg = V^F(r^, 3, Y) (10b) 

vAiere V^, 3, and Y depend oi wind speed. 

Alternatively, cne can use log-normal distributions to represent both 
conponents. In a recent study of atmospheric electrical conductivity this 
method (vhich included a log-normal cloud-droplet component, V^) \sas found 
to be reasonable . The NFS log-normal aerosol nunber density model is as 
follows : 

N(r) = N(r/2TT<^jj^) exp[-(J^n(r/rj^) )^/(20j^^)] (11a) 

cfj = ^n(cJ) (11b) 

v»here r^ = rjyg(S) is the dry (S = 0) radius 

g(S) = 0.8 (1-Ml-S)) , S > 0.4 (1^) 

g(S) = 1 (12b) 

The model parameters are given in Thble 1. 

Thble 1 . Log-normal aerosol parameters for the NFS model . Wind 

speed is u(m/s) and cloud liquid vater density is W(gm/]<g) . 
(Note that W depends cn height within cloud . ) 





r , lYn 
mo 


a 


N 


Cbnt . 


.027 


2.2 


10^ - 10^ 


Sea Salt 


.088 + .026 u 


2.2 


0.0012 


Cloud 


4 + 10 W 


1.4 


50 - 200 



A typical Junge value (A = 1 at S = .8) corresponds to N^ = 650 par- 
ticles per on^. A typical wind speed of 7 m/sec gives 55 sea-salt par- 

3 

tides per on with a nunber density mode radius of 0.27 Um. 

2.3 Relaxation of the Vfell-Mixed Assunption 

Considerable evidence exists to suggest that the veil-mixed assunption 
is an approximation. It is believed that this is due to the balance of 
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surface flux, entrainment and the boundary layer mixing time. R>r the 

qvsntity, X, let be the surface flux and be the flux at the top of 

the boundary layer (fluxes defined positive iJ^JVvard) . X can be q or 

aerosols. Define three velocities, W , W and W . 

e g o 



w = 

e 


-v*" 


(Ife) 


W = 
o 




(13b) 


W = 


raax(W , W ) 


(13c) 


g 


o e 




The ^ quantities are the junps in 


X at the two interfaces 




A 

o 




(14a) 


A = 
h 


^f -^mh 


(14b) 



vhere X^ is the sea-surfe.ce value, the value of X at Z = 0, the 



value of X within the mixed layer at Z = h and X^ the free atmosphere value 
of X at Z = h (Fig. 4) . The value of X within the mixed la^er at height Z 
is given by 

X(Z) = + ?(Wg/2Wg)^'^^ + (l-C)(W^2Wg)^'^^ (15) 

v^here C = Z/h, X^ is the average ''Viell-mixed" value and K is an integer (K 

= 3 is presently used) . Note that W = C u vhere C is the drag 

O u 

-3 -3 

coefficient. A typical value for W is 3 x 10 vshile W is 9 x 10 at 
u = 7 m/s. 

The application of Byi (15) to 9^ and q leads to a different (more 
gradual) height dependence of S in the mixed layer. The application to 
aerosols leads to a slight decrease in the total nunber density with in- 
creasing height. Both of these fectors result in a more gradual increase 
of extinction with increasing height in the mixed la^r. Therefore, f will 
be closer to 1.0. 



12 




X 



Figure 4. Flepresentation of the gradient of a meteorological variable 

in the boundary layer v^hen the \«ell-mixed assumption is relaxed. 
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2.4 Mixed-La >er Eepth 

The application of Bgn (6) is critically dependent i:pon a reliable 
estimate of h. In mid- latitudes doninated by high pressure (e.g., U. S. 

Vfest ODast) h is typically 200 to 800 m, in tropical regions h is typically 
1 to 2.5 hn. Note tlat h can vary a factor of two over a 24-hour period at 
a given location. It is hoped that h can be determined fron climatology 
(cortoined with synoptic scale classification) or ranote sensing of nearby 
cloud top temperatures. 

2.5 Synoptic CLassificaticn 

The method used for synoptic classification vas taten fron the Itefrac- 
tive Effects Guidebook (REG, CCMTEiIRFLT The. Man. 230-2-75) . Figure 5 is a 
sample synoptic chart frem the REG. The REG classes are given below in 
Thble 2. 



Thble 2. 


REG synoptic classes, typical duct heights, subsidence 
and air mass type (ST = si±>tropical , P = polar, T = 
tropical, M = modified) . 


TYPE 


h,m 


SUBSIDENCE 


A 


200 


-H- 


ST 


B 


600 


+ 


ST 


C 


200 


++ 


P 


D 


700 


+ 


MP 


E 


1700 


0 


T 


F 


3000 


0 


MT 


G 


4800 


- 


MT 


H 


100, 1700 


0 


T 


I 


High 


- 


T 


J 


2400 


0 


MT 


K 


no inversions 


0 


Ttemp 
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I 




Figure 5. Sample synoptic weather chart for the North pacific in Spring. 

The capital letters designate regions of different synoptic 
class 
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It must be remembered that the BEG designations of duct height are 
based primarily on discontinuities in the vvater vapor profile. The lowest 
water vapor disccntinuity may or may not correspond to the top of the mixed 
la^er. Rnr example, the E category data fron MACAT had insignificant water 
vapor junps at the top of the mixed layer (topically h ~ 250 m) . The major 
discontinuity in water vapor occurred above 1500 m, consistent wo.th the 
value given in Tbble 2. 

3.0 MAGAT Experiment 

The ^ferine Aerosol (feneration and Transport experiment (MACSAT) was 
held in Monterey in ^fely, 1980. The instrunented platform participating 
were the R/V ACANIA (NFS) and the Airborne Itesearch Associates (ARA) - 
Bellanca aircraft, iferosol sampling devices were installed on both plat- 
forms { Ihirall , 1981 ) . 

3.1 REG Q.assificaticns 

Ehch day during MACAT (4/29 - 5/09) was classified using the surfece 
pressure chart analysis from FNCX:. The results are given in Tbble 3. 



Thble 3. REG classes during MAC3AT. 



DATA 


CLASS 




4/29 


E 


(jEtont 


2300) 


4/30 


D 






5/01 


E 


(Front 


2300) 


5/02 


D 






5/03 


D 






5/04 


D 






5/05 


E 


(Front 0700) 


5/06 


D 






5/07 


D 






5/08 


E 


(Front 


1500) 


5/09 


E 
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3. 2 Aerosol Extinction Profiles 



Typical profiles for MAGAT D and E type classes are ^ovn in Fig. 6. 
The MAGAT profiles were divided into cloud- free and cloudy profile subsets. 
A profile ve.s designated as cloud-free if no cloud droplets vere encoun- 
tered during the profile measuranents . This does not imply that no clouds 
vere visible anyvhere but simply means that this particular sounding did 
not pass through a cloud (T^ble 4) . 

The average value of f for the cloud- free cases vas 0.90 with a stan- 
dard deviation of 0.22 and an uncertainty in the mean estinete of + 0.06. 
There vas no significant difference in the average f for D and E classes 
but the variance for D (based cn only four points) vas much less. Note 
that the optical depth la ' in Tbble 4 is the total actually measured by the 
aircraft (typically to a maximum altitude of 1.5 kn) , not the total atmo- 
spheric op^tical depth. The average mixed la^^r height under cloud- free 
conditions vas h = 0.29 kn vhile under cloudy ccxditicns h = 0.49 kn. This 
simply implies that the average lifting ccxidensation level (LCL) for the 
experiment vas no greater than 0.3 kn. 

3 . 3 Tropospheric Extinction 

A composite of all tropospheric extinction profiles frcm MAQVT (Z > h) 

is givQi in Fig. 7. Since the profiles do not indicate any significant 

height dependence cn this scale ( .4 to 1.5 kn) , the data vere conbined in a 

— —2 

number distributicn graph (Fig. 8). for D classes - 1.2 x 10 and for 
- —2 

E classes a_j_ = 3.0 x 10 • The standard deviation for LGT (a_j_) vas 0.5. 

The difference between for the two classes is probably not due to a 
difference in aerosol content but a reflection of the greater relative 
humidity (RH > 40%) present above the mixed layer for REG class E. 
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Figure 6a. terosol extinction coefficient (x = 0*49 micron) for a D class 
day 
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Figure 6b. Aerosol extinction coefficient (x = 0-49 micron) for an E class 



day 
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LGT (EXT(1/km)) 



Figure 7a. Tropospheric extinction (x = 0-49 micron) for Z > h and REG 



class D) 
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Height (km) 




Figure 7b. Tropospheric extinction (x = 0-49 micron) for Z > h and REG 
class E) 
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TS-ble 4. Data sumary for the 26 ladder profiles fron MAGIT. 

]i ' is the total optical depth to the maximun altitude 
m^sured and is the mixed la^er optical depth. 



CLEAR SKY 
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a >an ^ 
o — 
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— m- 


h,km 




REG 


4/29 


1 


.57 


.33 


.19 


.31 


.92 


E 


4/29 


2 


.46 


.40 


.37 


.55 


.68 


E 


4/30 


4 


.47 


.20 


.19 


.34 


.84 


D 


5/1 


5 


.14 


.074 


.018 


.16 


1.17 


E 


5/1 


6 


.23 


.11 


.082 


.24 


.67 


E 


5/4 


14 


.18 


.065 


.055 


.25 


.82 


D 


5/5 


17 


.17 


.040 


.019 


.15 


1.34 


E 


5/5 


18 


.13 


.035 


.012 


.15 


1.02 


E 


5/5 


20 


.10 


.029 


.020 


.18 


.9 


E 


5/6 


22 


.28 


.18 


.15 


.53 


.98 


D 


5/7 


23 


.46 


.14 


.12 


.23 


.88 


D 


5/9 


26 


.21 


.15 


.088 


.23 


.55 


E 




Average = 


.28 


.15 


.11 


D .34 


.9 














E .25 
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4/30 
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.5 






.45 




D 


5/2 


7 


.18 






.40 




D 


5/2 


8 


.15 






.44 




D 


5/2 


9 


.19 






.42 




D 


5/3 


10 


.21 






.52 




D 


5/3 


11 


.26 






.55 




D 


5/3 


12 


.19 






- 




D 


5/3 


13 


.56 






.50 




D 


5/4 


15 


.27 






.60 




D 


5/4 


16 


.16 






.42 




D 


5/5 


19 


.13 






.32 




E 


5/6 


21 


.30 






.60 




D 


5/7 


24 


.35 






.50 




D 


5/8 


25 


.20 






.90 




E 
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Ihe estimation of the tropospheric contribution, requires a scale 
length L vhich is either the exponential decay length or the height of the 
tropospheric aerosol ccraponent h^. If Vve tate a typical value of L = 7.5 
hn, then = 0*75 for class D and = .225 for class E for MACAT. Assum- 
ing a climatological mean for = .01 + .005 (Toon and Pollack, 1976) then 
.11/. 2 = .55 for D and = .11/. 35 = .32 for E cases. Alterna- 

tively, if ve assune h^ = 2 km (based cn REG data and sane aircraft ctoser- 
vaticns) then y^/u = .11/. 14 = .79 for Band u^/p = .11/. 18 = .61 for E 
cases. Thus, the mixed layer typically contributes 60% to 80% of the total 
optical depth for MAGAT. 

3.4 ^feriability Analysis 

The MAQ^T data have yielded sane estimates of the variability of the 
ccmpcnents in Bgi (6) v\hich can use for a simple error analysis 
(■Ihble 5) . 



Tcible 5. ^^riability contribution of the items of Hm (6) for MAG^T. 



E&rameter, X 


f 




^t 


h,kn 


Average , X 


0.9 


0.01 


0.09 


0.3 


4X 


0.2 


0.005 


0.07 


0.15 


ccntri. 

o 


Afp /h 
m 


fAPg/h 


f4p^/h 


fp Ah/h^ 
m 


Vhlue Aa 

o 


0.07 


0.015 


0.21 


0.17 


* 

This value of is the 

the aircraft. 


estimated 


total 


not the 


portion measured by 


A sun of the Vcuriances ( 


a^) of the 


four terms yields 


a total standard 


deviation a = 0.28 


for Oq (therefore 4a^a^ = 0.85) . 


3.5 Algorithm Evaliation 








Ihe MACAT data can 


be used to 


test Bgn 


(6) in a 


form optimized to the 



data set by using MACAT mean parameters for the tenns . If >e assune a_j_ = 
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0.01, H = 0.35 for class D and = 0.03 and R = 0.25 for class E, then 

a^' = .9(u’- 0.013)/1.35 for class D (16a) 

Oq' = .9(u'- 0.033)/l.25 for class E (16b) 

In lable 6 vve show a conpariscn of a^' with the measured for the 12 
MAGAT cloud- free profiles. 



T^ble 6. Evalvation of Etpi (6) as optimi 2 sd (Bcyi (16)) for the 
mean conditicxis during MAGAT. 



# 


Class 


u 

m 






a^, kn“^ 


“'/«o 


1 


E 


.33 


.29 


1.04 


.57 


1.82 


2 


E 


.44 


.40 


1.44 


.46 


3.13 


3 


D 


.20 


.19 


.49 


.47 


1.24 


4 


E 


.074 


.036 


.13 


.14 


.92 


5 


E 


.11 


.072 


.26 


.23 


1.13 


6 


D 


.265 


.052 


.13 


.18 


.72 


7 


E 


.04 


0 


0 


.17 


0 


8 


E 


.035 


0 


0 


.13 


0 


9 


E 


.029 


0 


0 


.10 


0 


10 


D 


.18 


.17 


.44 


.28 


1.57 


11 


D 


.14 
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Pn analysis of the 


data in 


Thble 6 yields /a^) 


= 1 . 07 and 





v»hich is consistent with the estimates of the previous section. 



4. 0 NPS Model Simulaticn 

The NPS aerosol model, as described in fection 2.2, has not been rua 

cn the MAGAT data sets. A preliminary simulation vas done for a typical 

climatological mean for open ocean conditions. The following parameters 

—3 

v\ere selected: N^ = 650 cm ,u=7m/s, h=.35 kn, RH(z = 10 m) = 80%, 

Tb = 15*^0, T = 13°C, Aq = -5, A0^ = 10 and = 0.1 on/s. This case 
yielded = 0.33 and f = 0.72 (Fig. 9) . The value of vas based on open 
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Figure 9. Sample NPS aerosol model simulation for typical open ocean 
climatological conditions (see text, section 3.4) 
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ocean measurements fron CEWCCM-78, JASIN and STOEX. Note that = 650 
cra”^ yields a_^ = 3 x lo”^/ vvhich is ccnsiderably lov^r than the values 
fbind during MAGAT. The R/V ACANIA measurements of the ccxitinental can- 

_3 

pcxient imply = 8 x 10 vdnch is in agreanent with the aircraft data. 

The model gave f = 0.77 viiich is ccnsiderably less than the value of 0.9 
found for the MAGAT data. 

5. 0 Conclusions 

5. 1 concept 

• Ihe optical depth is considered to be made ip of three atmospheric 

regions: mixed layer, tropospheric and stratospheric. 

• The relationship of surfoce visibility (Bcyi 6) is quite simple in 
the three layer context. 

5.2 MAGAT Evaluation 

• The REG synoptic categories vvere clearly correlated with certain 
features of the extinction profile (e .g . a_^_ and h) . 

• The value of f of 0.9 indioates that it is a smaller correction 
than anticipated from model simulation. 

• The tropospheric component of the optical depth varied about a 
factor of three . 

• Using parameters optimized to the data, Bgn (6) could predicot to 
roughly a factor of 2. 

5.3 NFS Aerosol Model 

• The model reproduced the general mixed layer and tropospheric 
extinction profiles . 

• Although the model allows some relaxation of the well-mixed 
assumption, it appears that the mixed layer gradients are slightly 
greater thain given by the model. 
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• The primary utility of the model in the future for this project 
may be in its ability to indicate the presence of a cloud- free (h < 
LCL) boixidary la^er. 

5.4 General Connents 

• The accuracy of the model (Bcyi (6)) is primarily limited by the 
variability of h and (see Thble 5) . 

• GLcbal application of Byi (8) using climatology only will probably 
allow estimates of 01 ^ to about a factor of 3. (Note ttet it is the 
accuracy of the Vfells-Munn-Fhtz model based cn in-situ wind speed 
and hunidity measuranents . ) 
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